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We demonstrate a method to determine dipole matrix elements by comparing measurements of 
dispersive and absorptive light ion interactions. We measure the matrix element pertaining to the 
Ca II H line, i.e. the 4 2 Si/ 2 «-)• 4 2 Piy 2 transition of 40 Ca + , for which we find the value 2.8928(43) eao. 
Moreover, the method allows us to deduce the lifetime of the state to be 6.904(26) ns, which 

is in agreement with predictions from recent theoretical calculations and resolves a longstanding 
discrepancy between calculated values and experimental results. 

PACS numbers: 37.10.Ty, 37.10.-x, 32.80.Qk, 03.67.Lx 


Methods for trapping and cooling single or few atoms, 
molecules or ions and manipulating them at the quantum 
level have opened up new avenues for precision laser spec¬ 
troscopy. In particular, quantum logic techniques [lj have 
enabled a new accuracy regime of timekeeping with opti¬ 
cal atomic clocks |2 . In contrast to atomic transition fre¬ 
quencies, dipole matrix elements and radiative lifetimes 
are still notoriously hard to determine at high accuracy, 
but are important for the quantification of black body 
radiation shifts of atomic clocks [3j, interpretation of as- 
trophysical spectra BUS], novel approaches for the search 
for physics beyond standard model 00 and for testing 
the accuracy of atomic structure calculations [5j. 

Regarding measurements of radiative lifetimes and 
transition matrix elements, established methods e.g. 
based on ion beams have been successfully complemented 
by novel techniques based on trapped particles. For 87 Rb, 
dipole matrix elements have been determined on the 10~ 3 
uncertainty level by diffraction in a condensate [§), while 
for the 6p 2 P°/ 2 state of m Yb + , the radiative lifetime 
has been measured by time-resolved counting of photons 
emitted from a single trapped ion m- A related tech¬ 
nique was used for neutral 1 ' 1 Yb in an optical lattice El. 

The species Ca + is widely used in quantum optics ex¬ 
periments, and its II H line led to the discovery of the 
interstellar medium [12]. For 40 Ca + , branching ratios 
between different decay channels have been determined 
at uncertainties approaching the 10 -5 level mm, and 
lifetimes of metastable states have been accurately mea¬ 
sured [T5j . The radiative lifetime of the 4 2 P 1; / 2 excited 
state of 40 Ca + has been determined to be 7.098(20) ns 
m by fluorescence measurements on a fast ion beam. 
However, this value disagrees with the most recent the¬ 
oretical calculation [5] by more than 11 standard devia¬ 
tions, while similar calculations for alkaline-like species 
from Li to Fr, Mg + , Ba + and Sr + are in good agreement 
with experimental results (see [8] and references therein). 



FIG. 1. (color online) Measurement scheme: (a) a laser 
off-resonantly couples two Zeeman ground state levels to an 
excited state with decay rate 7. (b) The coupling gives rise to 
a differential ac Stark shift Ag between the Zeeman levels and 
incoherent redistribution of population at rates R±. (c) Rel¬ 
evant energy levels in 40 Ca + . The shelving to the metastable 
3 2 D 5/2 state as well the 3 2 D 3 / 2 state acting as a population 
sink are indicated. State populations, which result in a bright 
(dark) detection event are marked by an open (closed) circle. 

In this work, we determine the radiative decay rates 
and the lifetime of the 4 2 P;i/ 2 excited state of 40 Ca + to¬ 
gether with the dipole matrix element of the 4 2 S;iy 2 f-t 
4 2 Pi /2 transition. The cornerstone of our scheme is the 
comparison between the dispersive and absorptive inter¬ 
actions, which occur upon driving this transition with an 
off-resonant laser, see Fig. [T] As the method is based 
on the discrete discrimination of atomic states of a single 
trapped particle mm, it is robust against many sys¬ 
tematic error sources which affect other existing methods. 

The off-resonant laser is characterized by its detuning 
A from the 4 2 Si / 2 transition, the Rabi fre¬ 

quency Ul and the relative amplitudes e q , which charac- 















2 


terize the circular (q = ± 1 , denoted ’±’ henceforth) and 
7 r {q = 0 ) polarization components. 

The dispersive interaction with the laser field, detuned 
by a frequency A from resonance, causes ac Stark shifts 
m [20] of the energy levels. Specifically, we are in¬ 
terested in the differential ac Stark shift A g between 
the two Zeeman-sublevels of the electronic ground state 
l'S'i/ 2 i m S = ±1/2), denoted henceforth as |f) and ||). It 
reads 


As 


I O 2 

- — ( f 2 - f 2 i 
3 4A “ + 


( 1 ) 


and is obtained by a spin echo measurement technique 

m- 

The absorptive interaction mediated by the same laser 
field manifests itself through inelastic Raman scattering, 
i.e. spin flips [ 221123 ] - We denote the spin flip rate from 
||) to |t) with R + and the rate in the inverse direction 
respectively with . For an optical field with arbitrary 
polarization they read [ 20 ] 


R± = IPS 


e± + £q f i 2 
9 4A 2 ’ 


( 2 ) 


where jps is the radiative decay rate from the P 1/2 to the 
S 1/2 state. The spin flip rates R± are experimentally de¬ 
termined by monitoring the population in |f) and ||) dur¬ 
ing the interaction. Note that the 7r-polarized field com¬ 
ponent drives the two spin-flip directions equally strong, 
whereas each of the the circular components drives only 
one pathway, corresponding to optical pumping in the 
limit of purely circularly polarized field component. The 
occurrence of elastic (Rayleigh) scattering is taken into 
account by the prefactor. The detuning A is chosen such 
that the conditions |A| fl, 7 ps for Eqs. [I] and [ 2 ] to 
hold are met. 

By combining Eqs. [l] and [2] we obtain the decay rate 
7 ps- 

Ips = 3A—, ( 3 ) 

with SR = R- — R+. All quantities on the rhs of [3] are 
experimentally accessible. 

As both interactions are driven with the same optical 
field, the Rabi frequency Q and the polarization ampli¬ 
tudes e q cancel out in Eq. [3] and notorious error sources 
are eliminated. These quantities are given by the elec¬ 
tric field amplitude of the off-resonant laser at the posi¬ 
tion of the ion, which reads E = X) 9 =± 1 |E|eqe g + c.c.. 
Here, e q are the spherical basis vectors. The polariza¬ 
tion amplitude e q define the coupling between Zeeman 
sublevels with mp — ms = q with effective Rabi fre¬ 
quencies scaled with the respective Wigner 3j-symbol, 
He q (-rrfs \ mp ) ■ Here H = 2|E| V/h is defined as the 
base Rabi frequency, with the reduced dipole matrix ele¬ 
ment V = (<Si/ 2 ||d||-Pi /2 )- 

We take an additional decay channel into account, for 
the case of 40 Ca + decay from the 4 2 P;w 2 to the 3 2 D 3 / 2 


state occurs at rate 7 pp. For that, we assume that only 
circular components are present in the beam, i.e. £q « 0. 
The rates at which population from |f) (||)) is sunk in the 
metastable 3 2 D 3 / 2 state are then proportional to R±: 

R^(1)d = 3(7pp/7ps)R_ (+ ) = 6 R_(+), (4) 

with the leak factor b , which is also accessible from the 
measured time-dependent populations. Together with 
7 pd, we obtain also the 4 2 P !/ 2 state’s lifetime r. Fi¬ 
nally, using 7 ps and the resonance wavelength A ps of 
the 4 2 Si / 2 ^I 2 Ri /2 transition, we can deduce the re¬ 
duced matrix element T>: 

X> 2 =2 T p S |±A> Ps . (5) 

We perform measurements on a single 40 Ca + ion stored 
in a segmented Paul trap [ 21 ] • The relevant energy levels 
and transitions are depicted in Fig. [l] A magnetic field 
splits the Zeeman sublevels of the electronic ground state 
|t) and It) by 27rxl3.7 MHz. The ion is Doppler cooled 
on the 4 2 S !/ 2 e->-4 2 P 1 / 2 (cycling) transition near 397 nm 
and afterwards prepared in either |f) or ||) by optical 
pumping. Now, the ion is illuminated with light near 
397 nm, detuned by A from the 4 2 Si / 2 f>4 2 Pp 2 transi¬ 
tion to induce both the dispersive and absorptive inter¬ 
actions, which allow for determining Ag and R±. Spin 
read-out is accomplished by shelving population from the 
|t) level to the metastable 3 2 D 5 / 2 state by means of rapid 
adiabatic passage (RAP) pulses [IBs 25]. This allows for 
discrimination between |t) and ||) as the linewidth of the 
4 2 <Si / 2 3 2 D 5 /2 quadrupole transition and the band¬ 

width of the RAP pulses are much smaller than the Zee- 
man splitting. Then, the ion is illuminated by laser fields 
driving the cycling transition and the 3 2 D 3 / 2 < r->-4 2 P 1 / 2 
transition near 866 nm, such that fluorescence is detected 
on a photo-multiplier tube if the ion has not been shelved. 
Conversely, the probability to not detect fluorescence, i.e. 
for a dark event, corresponds to the probability that the 
ion has been shelved from |f) to the metastable state. 

The off-resonant laser light is provided by an ampli¬ 
fied and frequency-doubled diode laser system, stabilized 
in both wavelength (5A/A < 0.5 • 10 -3 ) and intensity 
{SI/I < 0.5 • 10~ 3 ). The beam is aligned along the quan¬ 
tizing magnetic field and predominantly a + polarized, 
such that i?+ >> R-, and e 2 « 0, justifying the approxi¬ 
mation for Eq. [4] 

The ac Stark shift A 5 is measured with a spin-echo 
sequence, where a 7 r /2 pulse on the stimulated Raman 
transition between |f) and ||) is followed by a 7 r pulse, 
and another 7 t /2 pulse concludes the sequence. The de¬ 
lay time between the pulses is constant. During the first 
delay, the ion is exposed to a square-pulse of variable du¬ 
ration from the off-resonant laser. We probe 250 different 
shift pulse times, each with 150 interrogations. The pulse 
durations are spaced by 120 ns, such that a signal with 
about 40 oscillation periods is obtained. The differential 
ac Stark shift is the frequency of this oscillatory dark 
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FIG. 2. (color online) Dark event fraction versus time of ex¬ 
posure to an off-resonant laser pulse for initializations in |'|") 
and |t). Data points indicated by blue squares (red circles) are 
measured with the laser detuned by 12.03 GHz (13.94 GHz) 
from resonance. Solid lines show the fits to the model Eq. 
0. The inset shows the dark event fraction versus pulse du¬ 
ration for a spin-echo experiment. The blue line connects the 
datapoints, the red line is the resulting fit. The oscillation fre¬ 
quency directly corresponds to the differential ac Stark shift 
A s- 

event fraction (see inset of Fig. [2]). 

For the measurement of the spin-flip rates R± , the ion 
is exposed to pulses of variable duration up to 1 ms of 
the off-resonant laser after preparation. This changes 
the spin populations as depicted in Fig. [2j The time 
gap between initialization and readout is kept constant 
irrespectively of the scatter pulse duration to avoid sys¬ 
tematic effects. 

The spin-flip dynamics can be described by the rate 
equations: 

p t = -R_ (1 + b) p t + R + p± 

Pi = ~R+(l + b) pi + R-Pf. ( 6 ) 

The solution of Eqs. [6] is 

p^\t) = Ae - 2 ® ^cosh(^) — (1 + 6)<Si?sinh(^)^ 

p^\t) = ^e~2 m R +s inh(Rt/2), (7) 

with pi ' corresponding to initialization in |f),p^(£ = 
0 ) = 1, and for initialization in 1 1) ,p±(t = 0) = 1. 
We use R = (1 + &)(!?_ + R + ) and R 2 = R' 2 — 4&(2 + 
b)R-R+. 

For each measurement, the dark event fraction is deter¬ 
mined by probing the ion 2500 times for 30 fixed scatter 
pulse durations. Note that while the ac-Stark shift is 
given by the frequency of an oscillatory signal, the spin- 
flip rates are given by time-constants of exponentially 
decaying signals. Thus, the latter measurement requires 
significantly more data to attain the same level of preci¬ 
sion as the former. We extract the values for the differ¬ 
ential scattering rate SR and the leak factor b by means 


detuning from resonance [GHz] 



FIG. 3. (color online) Determination of the four detunings 
chosen for our measurements. For each of the 65 acquired 
datasets, we plot the quantity sign(A)As/51? oc A versus the 
optical frequency as measured by a wavelength meter. The 
zero crossing of a linear fit, depicted by the solid red line, 
reveals the resonance frequency with a standard measurement 
uncertainty of 2-7TX21 MHz. The insets show details of the 
datasets measured at the different detunings with the abscissa 
magnified 166 times. The shaded backgrounds indicate their 
uncertainties along the frequency axis. The arrows on the 
top frequency scale indicate the values of the four different 
detunings at -13.94, -12.03, 11.52 and 13.42 GHz. 

of a Markov chain Monte Carlo parameter estimation, 
taking into account the binomial statistics of the shot 
noise in the spin readout. State preparation and mea¬ 
surement errors are taken into account by modeling the 
measured dark event fractions as a linear transformation 
of the values p^\t),p!^\t) from Eqs. [ 7 J 

We repeated ac Stark shift measurements and spin flip 
rate measurements in an interleaved fashion in order to 
capture drifts of the laser intensity and thus the Rabi fre¬ 
quency f2. One measurement run consists of one spin flip 
rate measurement for preparation both in |f) and |)_), pre¬ 
ceded and followed by one ac Stark shift measurement. In 
total, 65 of such measurement runs were performed, at a 
total acquisition time of about 50 h. Four different values 
of the detuning A were used to show that no systematic 
effects with respect to this parameter are present. To de¬ 
termine the value of these detunings, we read the optical 
frequency from a commercial wavelength meter (High Fi¬ 
nesse WSU 267) with better than 10 MHz precision. We 
perform a linear regression of the measured A s divided 
by SR versus the corresponding optical frequencies to ob¬ 
tain the resonance frequency, as shown in Fig. [3] Com¬ 
pared to direct fluorescence spectroscopy, this technique 
mitigates effects like power broadening, Zeeman splitting 
or micromotion induced broadening, which would make 
a determination of the resonance frequency less accurate. 

Finally, we combine the measured quantities A, As 
and SR, using Eq. |3j to obtain the desired value 'yps- 
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FIG. 4. Comparison of the measured lifetime with theoret¬ 
ical and experimental results for P 1/2 - level lifetimes from 
references [H][Tfij| 26 j [27][ 28 l[29l[3(71 [3T][32l[33 l [34] 

Together with b this also yields 7 pd = §7ps, and the 
radiative lifetime is given by 1 /r = 7 ps(1 + |). 

Averaging over the data sets and taking into account 
the corrections discussed below, we obtain the resulting 
values: 7 ps = 27 r x 21.57(8) MHz and a —>-4 2 S !/ 2 

branching fraction of l/(l+|) = 0.93572(25). This 
leads to a lifetime of r = 6.904(26) ns and a value 
of 7 pd = 27r x 1.482(8) MHz. Using Eq. [5| we ob¬ 
tain a value for the reduced dipole matrix element of 
V = 2.8928(43) ea 0 . 

Experimental imperfections and model approximations 
lead to uncertainties and corrections for the value of 7 pg. 
These are summarized in Table [Tj Beyond the approxi¬ 
mation in Eq. |4j a possible eo polarization component 
has been taken into account, and yields corrections and 
uncertainties on the 10 -3 level. The uncertainty of the 
detuning A is determined by the finite precision of the 
wavelength meter, and by the uncertainty of the deter¬ 
mination of the resonance frequency, see Fig. [3] The 
random measurement errors for SR, determined by the 
amount of acquired data, is another significant contribu¬ 
tion to the uncertainty budget. Residual resonant light 
close to 397 nm resulting from imperfect laser switch-off 
causes relative corrections and uncertainties in the 1 CU 4 
regime. Beyond the rate equation model Eq. [ 6 j we in¬ 
clude a correction for the finite lifetime of the 3 2 D 3 / 2 
state or its depletion by residual light near 866 nm. A s 
can be measured by orders of magnitude more precisely 
than SR and b, however the precision is limited by drift 
effects, quantified by monitoring the ac Stark shift over 
the entire data acquisition time. Excess micromotion of 
the ion also causes a small systematic uncertainty via 
frequency modulation of the off-resonant light. Further¬ 
more, we include corrections due to the presence of the 
4 2 P 3 / 2 state about 27rx 6.69 THz above the ^ 2 R \/2 state, 
and due to the power-broadened Lorentzian lineshape be¬ 


Effect 

Shift x lO -4 

Unc.xlO 4 

Residual eo polarized light 

3.3 

2.9 

Uncert. of resonance frequency 

- 

1.6 

Stat. uncertainty of 8 R 

- 

1.5 

Residual near-resonant light 

-0.4 

0.5 

Wavemeter precision 

- 

0.4 

D 3/2 depletion 

0.3 

0.2 

Uncertainty of As 

- 

0.2 

Influence of micro motion 

- 

0.1 

Influence of P 3/2 state 

0.3 

< 0.1 

Residual line-broadening effects 

0.1 

< 0.1 

Total 

+3.6 

3.7 


TABLE I. List of relative corrections and measurement uncer¬ 
tainties for 7 ps- The specified values result from averaging 
over the 65 sets of measurement data. They are added in 
quadrature to obtain the resulting final uncertainty. 

yond the assumption |A| ft, 7 p$- A detailed discus¬ 
sion of corrections and uncertainties is presented in the 
supplemental material jS£l|. 

Our value for r is compared to previously reported 
values in Fig. |4j We find that our value agrees with the 
latest theory predictions [ 8 ] , while it is in substantial dis¬ 
agreement with the most recent experimental result 116) . 
Our value for the branching fraction is in agreement with 
the results from recent measurements M- To our knowl¬ 
edge, there is no experimentally determined value of T> 
so far. The value reported in our work is in agreement 
with the calculation from ||[j. 

Furthermore, using the measured value for 2?, 
we infer the reduced matrix element pertaining to 
the 4 2 Si /2 f->- 4 2 P 3 /2 transition, which results to 

(Si/ 2 ||d||P 3 / 2 ) = \f2V = 4.091(6) eao- Both reduced 
dipole matrix elements enter in the calculation of the 
blackbody-radiation (BBR) shift of the 4 2 S !/ 2 +)3 2 D 3 / 2 
and 4 2 S]y 2 +)3 2 D 5 / 2 quadrupole transitions, which are 
widely used testbeds for high precision laser spectroscopy 
155) . The BBR shift is among the major contributions 
to the uncertainties of clock transition frequencies, for 
species such as 43 Ca + [35], Sr + [5?], Sr [58] [Mj and Yb 
[401 . These are used for state-of-the art optical clock 
experiments, and some of them are discussed as new SI 
frequency standards. Our result adds to existing work 
[4T1 42, 03] validating computational methods used to 
predict the BBR shift in optical clocks. 

In summary, we demonstrate a novel method for the 
measurement of dipole matrix elements, which works de¬ 
spite the presence of additional decay channels. We at¬ 
tain an uncertainty on the 10 ~ 3 level. All major error 
sources can be potentially mitigated, such that measure¬ 
ments of radiative decay rates at unprecedentedly low un¬ 
certainties in the 10~ 4 regime appear within reach with 
current technology. Our method is applicable to atom 
and ion species which allow for preparation and readout 
of Zeeman or hyperfine sublevels. 
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In this document, we present details on how to avoid, suppress and characterize corrections, systematic 
and statistical errors in our measurement scheme to determine dipole matrix elements. In Sec. 1, we 
give details on how data for the different quantities was acquired, and numerical values for the resulting 
quantities are listed. Sec. 2 explains details of the experimental setup and methods which are related to 
error mitigation. In Sec. 3, we discuss all sources of corrections and errors in detail and list numerical 
values for these. Note that we refer to expression x from the main manuscript as Eq. Mx. 
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1 Measurements to determine 77 pd and r 

According to Eq. M3, the natural linewidth 7 pg is determined by the detuning from resonance A, 
the differential ac Stark shift A 5 and the differential scattering rate SR. These quantities have been 
measured as follows: 

• Detuning A: 

The detuning Ap is determined by measuring the the absolute frequency f :j of the scattering light 
for each detuning j. The number of detunings D = 4 (two on the blue side and two on the red side 
of the resonance) and the resonance frequency f res of the S — P-transition, which is determined 
as described in section 3.2. 

• Scattering rates P+ and P_ 

We have acquired a total number of 65 scattering curve pairs, distributed among four detunings 
Aj, such that for each detuning A) = {18,17,14,16} measurements were done. Each pair consists 
of one curve for initialization in |f), and one in |{). For each of these data sets, 30 different pulse 
times when the off-resonant beam is switched on were probed, ranging between 0 ms and 1 ms. The 
pulse times are unevenly distributed along the time axis to provide a more dense sampling for short 
pulse times. The 30 pulse times are probed in chunks of 25 repetitions, in a random permutation 
to avoid systematic effects. This is repeated 100 times, such that in total 2500 interrogations per 
pulse time are performed. For zero pulse time, in total 5000 interrogations are performed, which 
allows for a better accuracy in the determination of the readout parameters, c.f. Sec. 2.3. The 
measurements for each initialization level are performed within one sequence. 

The polarization of the off-resonant light is adjusted to be predominantly cr + , i.e. P+ P_, 
and it is not changed throughout all measurements. The two curves of each curve pair are jointly 
fitted to the model function. The fit provides the values R+,ij and R-,i,j for curve pair (i,j). j 
identifies the detuning, and i identifies the number of measurement for that respective detuning. 

• Differential Stark shift A 5 

For each scattering curve pair (i. j). right before and right after the scattering measurement, Stark 
shift measurements are carried out, the mean of which, A s,i,ji is considered to be the Stark shift 
belonging to measurement (i, j). 

For each scattering curve pair, ”/ps,ij is calculated. We take the mean of those, and arrive at: 

7 p S = 2tt x 21.48 MHz (1) 
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Taking into account all corrections from section 3, that value reads: 


7 PS = 2tt x 21.57 MHz (2) 

For each measured scattering curves the branching parameter b, :i is also given by the same fit 

mentioned above. With the mean of all bij, b , the natural linewidths ~/ps and 'ypD are related as follows: 

7 pd = — 7 PS — 2-7T x 1.48 MHz (3) 

3 

The total lifetime of the Pi/ 2 -state r is finally: 

t = --- = 6.90 ns (4) 

7 PS + 7 pd 

The resulting values for these quantities are shown in table 1. The numbers are given for the four 
detunings A i; each averaged over the iV* respective datasets. 


i 

Ai/2n (GHz) 

Ri,+ {s~ l ) 

Ri,- (s -1 ) 

A s>I /2^ (MHz) 

h 

7PS,i/27r (MHz) 

1 

-13.94 

4474(78) 

47(7) 

1.3753 

0.2063(48) 

21.45(33) 

2 

-12.03 

6011(80) 

63(11) 

1.58819 

0.2080(53) 

21.52(30) 

3 

11.52 

6514(125) 

101(7) 

1.64198 

0.2056(45) 

21.49(39) 

4 

13.42 

4762(40) 

76(8) 

1.39819 

0.2067(52) 

21.48(17) 


2 Experimental considerations 

2.1 Post-selection 

The data is post-selected in such a way that chunks of 25 interrogations where the ion is not found 
sufficiently bright after reset - indicating an impairment of Doppler cooling - are discarded. Similarly, 
chunks where the shelving was impaired were discarded, which can be identified by spuriously low dark 
count fractions. In total, less than 1% of all raw datapoints have been removed by this post-selection 
scheme. 

2.2 Pulse area stability 

Our measurement scheme relies on a constant optical power and polarization from the off-resonant laser 
near 397 nm at the location of the trapped ion. Several measures are employed to guarantee the stability. 
The laser light derived from a doubled, amplified laser system (Toptica TA-SHGpro) is coupled into a 
single mode fiber. At the fiber output, a polarizer filters undesired polarization components. After the 
polarizer, a small fraction of the optical power is picked and measured on a photodetector. The signal 
is used by a commercial intensity stabilization unit (NoiseEater, TEM Messtechnik), which acts on an 
electro-optical modulator before the fiber input. This way, intensity and polarization after the fiber 
output are actively stabilized. For switching, we employ an acousto-optical modulator (AOM) 1 , which 
exhibits no measureable 2 beam pointing drift when switching between no power and maximum rf-power 

1 I-M110-3C10BB-3-GH27, Gooch&Housego 

2 Measured with a position sensitive detector PDP90A, Thorlabs 
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(2 W). The rf amplifier supplying the AOM is mounted on a temperature-stabilized water-cooled plate 
to avoid temperature-change-induced power drifts. The laser beam is focused such that zero and higher 
diffraction orders coincide at the position of the trapped ion, i.e. residual AOM induced beam pointing 
fluctuations are mitigated. Further mitigation is ensured by a relatively large beam waist of about 
30 /im, and a careful positioning of the ion at the waist center by maximizing the ac Stark shift. 
Diffraction angle and efficiency of an AOM weakly depends its temperature, and thus on the time- 
averaged rf-power. This is kept constant by fixing all pulse sequences for measurement of the ac Stark 
shifts and scattering rates to the same duration, and inserting additional rf pulses after readout to keep 
the average switched-on time constant. 

2.3 State preparation and readout (SPAM) errors 

Taking into account imperfect initialization, the experimental signals for the scattering rate measure¬ 
ments are modeled by a linear transformation of the populations from Eq. M6: 

4 T) (*) = 0- + (! - e t) (*) + 0- - a t)Pt J ’ ) (*)) 

4 40 (*) = o- + c 1 - e t) ( a t4 4) w + - «t)4 T) w) > ( 5 ) 

where the < 1 and cq. > 0 characterize the state preparation (pumping) fidelity, e| (e-f) is the error 
probability that the ion is spuriously detected as being in state |f) (l-D), when it actually is in state |f) 

(It))- 

3 Uncertainties and corrections of the measured variables 

We carefully distinguish statistical and systematic errors, and systematic shifts. A statistical error of 
a given quantity is denoted by e(-), while a systematic error is denoted by er(-). A systematic shift is 
denoted by S(-), and it is accompanied by a systematic error. All shifts and uncertainties are reported as 
relative quantities. Some error sources are found to be small as compared to the dominating ones, in this 
case an estimation of the magnitude is given, however these error sources are not quantitatively taken 
into account in the total error. If this is the case, it is explicitly stated in the corresponding section. 

As our total error is on the 10 -3 level, we do not take corrections and errors below the 10~ 4 level into 

account in the final error budget. In the cases where we report mean values e, cr, these quantities result 
from averaging over all data sets. In these cases, the individual values for each data-set rather than the 
indicated mean values are used for the computation of the final values. 

3.1 Statistical Uncertainty of SR and b 

For each curve i, R±,i.j andfr are determined by a fit of the depolarization data to the model Eqs. 5. 
The uncertainties of that fit are extracted by marginalization of the data generated by the used Markov 
Chain Monte Carlo algorithm. It is verified by parametric bootstrapping: For fixed R± and b values, 
30 instances of simulated measurement data with binomial noise are generated and fitted in the same 
way as the original data, each fit yielding result values for R± and b and their statistical errors. The 
standard errors determined from the actual data are in agreement with the standard deviation of the 
mean values found for the 30 bootstrapping datasets. Note that cross-correlations between uncertainties 
of R± and b and the additional parameters introduced for characterizing SPAM errors, c.f. Sec. 2.3, are 
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fully taken into account due to the marginalization of the Monte Carlo data. 

The individual standard error of SR for each dataset i, j is used to determine the total uncertainty of 
7 P 5 , see section 4. The individual standard errors of b affect the total uncertainty T pjj and t, c.f. Eqs. 
3 and 4. The means over all datasets of the relative statistical errors are 


e(5R) = 12.1 • 1CT 3 

e(b) = 24.8 • 1CT 3 (6) 




Figure 1: Comparison of expected and actual spreads of measured dark count fractions: We plot the 
expected population standard deviation Eq. 7 (solid black line), computed using the total estimated 
dark count fractions p)j’ (t) and N = 25. The red points are the sample standard deviations of the 
dark count fractions from the 100 different acquired chunks of N =25 interrogations. The solid grey 
lines show 2 a intervals for the expected standard sample deviations. A significant fraction of outliers 
would indicate additional drifts and fluctuations. The data for p^" (t) is shown in the left panel, the 
data for pi ; (f) in the right panel. 


The validity of the statistical errors relies on the assumption that their dominant source is the readout 
shot noise of the spin measurements. Possible drifts or fluctuations of the optical power or polarization 
during a depolarization rate measurement would corrupt this assumption. We verified that no such 
effects are present on a significant level in the following way: For a given data-set, for each exposure 
time t, the measurement data for pi (t) was taken in chunks comprised of N =25 spin interrogations. 

For each chunk, coarse estimates for p?’"^(f) are obtained by dividing the number of dark count events 
by the number of interrogations N. The sample standard deviation of these values over the set of chunks 
is compared to the expected population standard deviation for binomial statistics, 

e (pj t4) (f)) = (l-pf’ W (t)))/W (7) 

The results are in agreement with the assumption of binomial statistics, we can therefore exclude 
additional drifts and fluctuations at a significant level. An example is shown in Fig. 1. 
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3.2 Detuning A 

The detunings read A j = fj — fresi with f res being the so far unknown resonance frequency of the 
transition. The precision of the detunings A ? are given by the measurement precision of the absolute 
laser frequencies fj of the four detuning values. These are measured by a commercial wavemeter 3 . 
The resulting uncertainties have been quantified by determining the standard deviation of a time series 
of wavelength measurements of a frequency-stabilized laser, which has a linewidth orders of magnitude 
smaller than that measured standard deviation. These uncertainties amount to: 


e wm (A 1 ) = 0.4-10- 3 
e wm { A 2 ) = 0.5 • 10- 3 
e W m{ A 3 ) = 0.6 • lO" 3 

G«m(A 4 ) = 0.4 ■ 10 -3 (8) 


To determine the resonance frequency f res , we note that the ratio of ac Stark shift and differential 
scattering rate is proportional to the detuning: 


Ho 


1 Asyj 


3 SR 


OC A ; 


'hJ 


(9) 


The minus sign holds for negative detunings. We perform a linear regression of this quantity with respect 
to the wavelength meter frequency fj from which the zero-crossing / res , he. the resonance frequency, is 
obtained. The accuracies of the detunings A j are limited by the accuracy of the linear regression used 
to determine the numerical value of fres- As f res is used to compute all detunings, it is a systematic 
error. The values are: 


Vfit(A-i) = 1.5 • 10 -3 
o’ / it(A 2 ) = 1.8 • 10 -3 
( r /lt (A 3 ) = 1.8-10- 3 

u /lt (A 4 ) = 1.6-lO" 3 (10) 

3.3 Residual near-resonant light at 397 nm 

Two laser sources provide spurious residual light weakly driving the cycling transition. 

• Doppler cooling laser: The laser used for Doppler cooling and fluorescence readout is switched 
by an acousto-optical modulator 4 in double-pass configuration. An estimate of the residual 
photon scattering rate due to imperfect switch-off is obtained by initializing the qubit in |f) or 
IJ.), either before or after a 10 ms wait time before readout. This yields spurious spin flip rates of 
aH 6 ) =1.1(6) Hz and R + es ' ) =0.7(4) Hz. These rates are taken into account for two of the acquired 

3 WSU, High Finesse GmbH 

4 I-M110-3C10BB-3-GH27, Gooch&Housego 
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datasets, by modifying the model Eqs.M 6 , incorporating the fact that these persist throughout 
the entire 1 ms interval. From fitting with the modified model, obtain the following corrections: 

S 397 (SR) = -4.1 • 10 " 4 

S 397 {b ) = -8.3-10 - 4 (11) 

and systematic errors 

<j 397 (5R) = 4.6 • 1CT 4 

a 397 {b) = 6.4 -lCT 4 (12) 


• Off-resonant laser: The off-resonant light near 397 nm is obtained from an amplified, frequency 
doubled diode laser system 5 . The amplified fundamental beam seeding the doubler has a broad 
background of amplified spontaneous emission (ASE), extending over about 30 nm, equivalent to 
Samp ~ 27t- 57 THz. The suppression with respect to the coherent light is specified to be \ =-40 dB. 
One might argue that this background leads to additional UV light near the atomic resonance, 
which would effectively shorten the measured lifetime. We give a pessimistic estimate of this effect, 
assuming only % =-30 dB ASE suppression, perfect phase-matching for sum-frequency generation 
(SFG) of ASE and laser mode photons, perfect AOM diffraction even for incoherent light and 
that a spurious component exactly hits the atomic resonance. With a SHG cavity linewidth of 
S c = 27T-500 kHz, we obtain a power ratio of the spurious to the main SHG component of 


P A 

rs P ur ~ ° c ~ q . in- 15 

PSHG X S amp 


(13) 


This strong suppression even under pessimistic assumptions leads us to the conclusion that this 
effect does not have to be taken into account. We still perform an experimental verification, where 
we probe depolarization after initialization in |1), at a fixed off-resonant pulse-time of 50 p,s, 
corresponding to about 20% depolarization. We scan the laser frequency over a range of 1.3 GHz, 
which is more than one free spectral range of the SHG cavity. The scan speed is about 2 MHz/s, 
such that one data point integrates over about 4 MHz, which is slow enough to resolve the atomic 
line. No resonant features were observed in the resulting depolarization signal. 


3.4 Residual 7r-polarized scattering light 

The k-vector of the scattering light is aligned parallel to the magnetic field, so ideally only cr-transitions 
can be driven. Alignment errors of the beam might be the cause for residual 7 r-polarized light, which is 
not accounted for in Eq. M4. Taking into account an the pumping to the D 3/2 caused by additional n 
polarized light, the rate equations Eqs. M 6 are extended, including an additional parameter Rd- 


Pf = -R- (1 + b) pf - R d pf + R+pi 
Pi = -i?+(l + b) Pi - R d pi + R-P\. 


with 


5 TA SHGpro, Toptica Photonics AG 


R D = b -^R 

2 el 


+ ■ 


(14) 

(15) 
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The solutions of the rate equations are of the same form as for Eq. M6, however with altered coefficients 


R = (l + b)(R- +R+) + R d (16) 

R 2 = R 2 -4(2 + b)(R D + bR_)R+. (17) 

If the data is now fitted using this modified model, we obtain values of about Rd =3(3) Hz, and altered 
values for SR, b. The mean value of the deviations of SR, b are used to compute the corrections: 

SJ5R.) = 3.3- 10” 3 

5^{b) = -2.4- 1(T 3 , (18) 

and the spread of the deviations within each data-set determines the systematic error 

a n (SR) = 2.9 • 1(T 3 

a n (b) = 2.6 • 10 -3 (19) 


3.5 Differential Stark Shift A g 



Figure 2: The red points depict the relative deviations of the ac Stark shift values A s,i,j from their mean 
for all four measurement sessions. The blue curves are the B-spline functions defined by the points. The 
RMS deviation of the points from the spline is used as the statistical error. 


We carried out the bootstrapping method from Sec. 3.1 for determining the uncertainty of the ac Stark 
Shift A S - The result is that the uncertainty from the binomial readout statistics of the measurement 
would lead to insignificant uncertainties on the 10 - ' level. Systematic effects, i.e. temporal drifts, are 
clearly dominant. 

A s i j is for each dataset i,j is determined by the mean of two measurements, taken before and after the 
scattering rate measurement yielding SRij. The drift of the ac Stark shift drift over the measurement 
time each detuning A,; is slow and therefore mostly resolved, thus the assumption of a linear drift between 
two consecutive stark shift measurements is justified. To account for residual drifts on a timescale which 
is not resolved, we use the measured Ag values as control points to define a B-spline function as shown 
in Fig. 2. The RMS of the difference of the measurement result and the corresponding B-spline value is 
treated as the overall inaccuracy e(A s,i,j) of all As values pertaining to detuning A*. The mean value 
of those relative uncertainties reads 

e(Ag) = 3.0 • 10" 4 (20) 
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3.6 Zeeman shift 


We do not take the Zeeman shifts of the atomic levels (A z — 27t- 13.34 MHz splitting between 
|t) and |t)) into account. However, as the atomic resonance frequency is determined via the ac 
Stark shift caused by the predominantly cr+ polarized beam, the detuning A specifically refers to the 
\S 1 / 2 ,mj = —1/2) \Pi/ 2 ,rnj = +1/2) transition including its Zeeman shift. Corrections are only 
necessary for the transitions driven by the weak er_ and 7 r polarization components. The corresponding 
systematic errors are estimated to be 


<tz(SR) 

cz(As) 



< 2 • 1(T 5 


— <3-10 " 5 

i?+ A 


( 21 ) 


3.7 Influence of Micromotion 


Both the scattering rates f?+, i?_ and the Stark shift As are altered in the presence of micromotion 
along the propagation direction of the off-resonant beam. The micromotion causes effective frequency 
modulation (FM), characterized by the trap drive frequency Q r f = 27T-25 MHz and the modulation 
index fj. We obtain 


^(mm) p 

e 

- 7 ps- 

A ( mm ) - 
^ S 

. 4 - 

Q 


e ± 


2 +°° 


E 


n 2 


n——oo 
2 +°° 


E J «^) 


4(A + nPl r f) 2 

n 2 


4(A + ntt r f) 


( 22 ) 

(23) 


By measuring the the shape of the resonance fluorescence line by recording fluorescence photons at low 
power and across resonance, an upper bound for the modulation Index of /? =3.5 is obtained. Note that 
the micromotion-induced frequency modulation leads to symmetric FM sidebands, therefore the effects 
cancel out to first order. Due to the A _ 1 (A~ 2 ) scaling behavior of the ac Stark shift (scattering rates), 
both quantities are slightly increased by FM components on the side towards the resonance, such that 
the total effect partially cancels out. Thus, we give an estimate for the systematic error of the resulting 
7 ps rather than for As and SR. Using realistic values of D 2 > e + > e - and A, we obtain 

ammhps) = ^ - SR) - ^ ( A 5 mm) - A s) » 1.0 • 10 ” 4 (24) 


3.8 Influence of near-resonance-effects 

The expressions for SR Eq. M2 and As Eq. Ml are only valid if |A| 7 ps,U. Our experimental 

conditions are 7ps/|A| < 2 • 10 ~ 3 ,D/|A| < 4 • 10~ 2 , we therefore investigate how the approximations 
affect the results. 

The measured , i.e. actual scattering rate SR( meas ^ is smaller than the far-off-resonance approximation 
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Eq. M2. Including power broadening and finite natural linewidth, the scattering rate is given by 


SR (meas) = ^!+_ 


n 2 


4A 2 + § ft 2 + 7 2 s 


IPS 


2 2 
e + “ e --n 2 


4A 2 


§» 2 +7ps 

16A 4 


= SR + SR' 


We thus have 


SR' = -6R 


I^+Yps 


4A 2 


The correction is obtained by invoking the expression Eq. M3 for the decay rate 'yps- 

Q A 

7 PS = — (SR^ mea ^ - SR') 

A s 


3A 

As 




such that 


We obtain 


Sbroad^SR) — 


SR' 


g R (meas) 


2A S /A 


(25) 

(26) 


(27) 

(28) 


SbroadiSRi) = 2.0 • 10 - 4 
Sbroad{5R 2 ) = 2.6 • 10 - 4 
Sbroad(SR 3 ) = 2.9 • 10- 4 

SbroadtSRi) = 2.1 ■ 10~ 4 . (29) 

For As, no expression corresponding to Eq. 25 is available, we therefore extract values from a numer¬ 
ical simulation using a three-level generalized optical Bloch equation. For Sbroad{As) = (A^ meas ^ — 
As)/Yg leas \ we obtain 


Sbroadi&S, l) = 1.0'10- 4 

Sbroad{As } 2) = 1.4 • 10^ 4 
Sbroad{As t 3) = 1-4'10 _4 

Sbroad( As, 4 ) = 1.0 • 10- 4 . (30) 

3.9 Finite .D 3 / 2 -lifetime and residual 866 nm repump light 

In the model Eqs. M6, the lifetime of the -D 3 / 2 -state is considered to be infinite. The most recent 
experimental value is 1176(11) ms (Kreuter et. al., PRA 71, 2005). Residual light near 866 nm also 
depletes population accumulated in the D 3 / 2 state at a rate Rdep, effectively reducing its lifetime. 
We measure the residual power in the imperfectly switched off beam, using a fibre-coupled avalanche 
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photodiode, to about 6 pW. In the switched-on state, the optical power is 86 /xW. For observing the 
amount of resonance fluorescence versus optical power at 866 nm, we infer a saturation parameter of 
23(5). This leads to a saturation parameter in the switched-off state of S^P =1.6(8)T0 -6 . The 
depletion rate of the D 3 / 2 state is then 


0(866) 

R (dep) _ TPS ^off 

D3 /2 2 1 + 4^/tV 


(31) 


Here, S is the average detuning of the sub-transitions between the P\/ 2 and -D 3/2 manifolds, considering 
the Zeeman splitting (see Sec. 3.6) and the repump laser to be tuned onto the blue side of the entire 
manifold. We obtain a laser induced depletion rate of 1.2(1.0) Hz. Thus, the D 3 / 2 state is depleted at 
a total rate of about Rp 3 % =2.2(1.0) Hz. 

For estimating the impact of the depletion on the result for SR, we modify the rate equations Eqs. M 6 
to be 


Pf = -R-0- + b)p t + R+Pi + \R^(l - - Pt) 

Pi = -R+(l + b)pi + R-Pt+^Rpyl{l-Pi-Pi). 

Using these for example data in our fit procedure, we infer the systematic shifts 

S-yosMSB.) = 3.2-lCT 4 
SyD3/2(b) = 1.6-HT 3 

and systematic errors 

^m/2{SR ) = 1.6-10- 4 
^D3/2(b) = 7-l(T 4 


(32) 


(33) 


(34) 


3.10 Influence of the P 3 / 2 state 


The P 3/2 is separated from the P\/ 2 manifold by the fine-structure-splitting of A ps ~ 27r • 6.69 THz. 
The coupling strength H 2 is increased with respect to the Pi / 2 state by a factor of 2, given by the ration 
of the squared reduced dipole matrix elements, which we indicate explicitly. The presence of these levels 
affects both the Stark-slrift and the scattering rates. For the scattering, we approximate the decay rate 
to be the same as for the 2 P\/ 2 state, and we take only the predominant er+ polarization component 
into account, and we neglect decay into the D 3 / 2 and D 5 / 2 states. The additional scattering rate then 
reads: 


d (P3/2) 2 e\ 2 fl 2 1 2 VL 2 

R + R , 7PS -___^« 7PS __ 


(35) 


Note that the for the P 3 / 2 state, due to the different Clebsclr-Gordan coefficients, the branching ratio 
of decays with spin-flip to decay without spin-flip is reversed with respect to the P 1/2 state, however 
the scattering rate has the same prefactor as in Eq. M2. The excitation channel Si/ 2 ,mj = +1/2 -+ 
P 3 / 2 ,mj = +3/2 contributes only to elastic scattering. Interference effects between different pathways 
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do not play a role as the emitted photons are distinguishable either by wavelength or by polariza¬ 
tion/emission direction. The correction to the ac Stark shift reads: 


A 


(-P3/2) 

S 


V 2 6 ) 4(A fs - A) 


e 2 + 2 n 2 
3 4A F s 


(36) 


For the P 3 / 2 state, the transition = +1/2 —► P 3 / 2 ,mj = +3/2 contributes to the differential 

ac Stark shift and has a larger Clebsch-Gordan coefficient than the transition S'i/ 2 ,toj = —1/2 -+ 
P 3 / 2 ,toj = +1/2. The reversed sign has to be accounted for as we only measure absolute values of As- 
For A > 0, the magnitude of A g is increased, for A < 0 it is decreased. 

We finally obtain 


<5p3/2(As,i,j) 


Aj 

3A fs 


1.0- 1CT 3 


Aj 

2n ■ 10 GHz 


(37) 


fip3/2(3Rs,i,j) ~ 

Note that the Sp 3 / 2 {As } ij) depend on the 
small magnitudes of these corrections. 


A/ 

3A fs 


1.5 • 10 


-6 


A,. 


(38) 


2n ■ 10 GHz, 

sign of Aj. We do not include systematic errors due to the 


3.11 Fit bias 

Parameter estimation from measurement data might be affected by a systematic bias, depending on the 
used estimator. We generated artificial measurement data from the model Eq. M6 with fixed values of 
R±,b with binomial errors. We find that for increasing number of measurements, i.e. vanishing noise, 
the fit results for R±,b converges to the preset values. We thus do not take any bias into account in our 
error budget. 

3.12 SPAM errors 

For the determination of the differential scattering rate (5P, SPAM errors are taken into account by em¬ 
ploying the model introduced in Sec. 2.3 for the parameter estimation. The uncertainty of the additional 
parameters and their cross-correlation with R±,b is already taken into account, as the statistical error 
e(SR) is obtained by marginalizing over the sample data from the Monte Carlo parameter estimation 
algorithm. 


4 Total uncertainties and corrections of 7 pp,7pi? and r 

The uncorrected values for jps for each detuning Aj are obtained from averaging: 
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The corrected values are obtained by adding all relative corrections: 

'yPS,i/'YpS,i = 1 T ^397(6-R) + Abroad.(8Ri) + SyD3/2(bR) + S n (6R) 

- Sbroad(As,i) - 5p 3 / 2 (As,i) 

h/K = 1 + fi397(b) + S^ D3 / 2 (b) + d„(b) ( 40 ) 

The statistical errors are obtained from error propagation: 

c(7 PS,i,j) 2 = e(A;) 2 + t{As,i,j) 2 + e (bRi,j) 2 

Ni 

€ (bi) 2 = ^fY^d b i,j) 2 ( 41 ) 

* 3 =1 

For taking the systematic errors into account, we face the problem that most of their correlations are not 
quantitatively available. However, it can be reasoned that they are due to well-distinguished physical 
effects and thus influence the resulting quantities in ways which are hardly related to each other. 6 7 . 
According to JCGM guidelines ' we therefore neglect possible correlations and add them in quadrature: 


O’(7 PS,i ) 2 = Pmmhps ) 2 + o/it(Aj ) 2 + a 397 (SR ) 2 + cr 7 D 3 / 2 (AR ) 2 + a n {SRi ) 2 

a (bi ) 2 = 0397 ( 6 )“ + cr lD 3 / 2 (b ) 2 + crx(bi ) 2 (42) 

The final values are obtained by a weighted average over the four results: 

7 ps = ^ AT y Nilps,i 

‘ m 

In Table I in the main manuscript, the mean values of all statistical and systematic uncertainties of -)ps 
are given. When added in quadrature, they yield the total uncertainty. In the following calculation, 
we pursue a more detailed approach, where the uncertainties of the 65 measurement runs are evaluated 
individually, before the mean is taken. The difference of the results however turns out to be insignificant 
( < 1 x 1CT 4 ). 

The total statistical errors are given by 

e(7ps) 2 = 1 2 Y, e bPS,ij ) 2 

ij 

e(6) = (E iNi) 2 ^ Nie{bi) ■ (44) 

6 Note that all systematic errors depend on A, except for afn(Ai) - the error of A itself. The resulting corrections from 
these correlations are on the 10 —6 level and are thus insignificant. 

7 Joint Committee for Guides in Metrology, Evaluation of measurement data - Guide to the expression of uncertainty 
in measurement , p. 62 (2008), http://www.bipm.org/en/publications/guides/gum.html. 
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For the total systematic errors, we assume complete correlation between experiments at different detun¬ 
ings and average them linearly: 

v(lPs) = y N N i a hPS,i) 

v(b) = y 1 ]sr ^2 N i<p( b i), (45) 

The total measurement uncertainties result from adding statistical and systematic errors in quadrature: 

£(7ps) 2 = e(jps) 2 + c(7ps) 2 

£(b ) 2 = e(b ) 2 + a(b ) 2 (46) 

For the radiative decay rate to the -D 3/2 state, 7 pp> = (5/3)7 ps, we thus find 

£(1pd) 2 = e{b ) 2 + c( 7 ps ) 2 + cr(b ) 2 + < 7(7 ps) 2 (47) 

and correspondingly for the total lifetime r = 1 /( 7 ps + 7pp): 

£(t ) 2 = e (7 ps) 2 + < 7(7 ps) 2 + + (e{b ) 2 + a{b) 2 ) . (48) 
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